A novel antiproton radial diagnostic based on octupole induced ballistic loss 
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We report results from a novel diagnostic that probes the outer radial profile of trapped antipro- 
ton clouds. The diagnostic allows us to determine the profile by monitoring the time-history of 
antiproton losses that occur as an octupole field in the antiproton confinement region is increased. 
We show several examples of how this diagnostic helps us to understand the radial dynamics of an- 
tiprotons in normal and nested Penning-Malmberg traps. Better understanding of these dynamics 
may aid current attempts to trap antihydrogen atoms. 



I. INTRODUCTION 

Cold antihydrogen atoms (H) were first produced 
by the ATHENA collaboration [1], and, shortly there- 
after, by ATRAP [2] at the CERN Antiproton Decel- 
erator (AD) [3] in 2002. They were produced by mix- 
ing positrons (e+) and antiprotons (p) held in Penning- 
Malmberg traps. Such traps use a solenoidal axial mag- 
netic field Bz to provide radial confinement, and elec- 
trostatic wells to provide axial confinement. Penning- 
Malmberg traps confine only charged particles and, con- 
sequently, do not confine neutral H atoms. 

The current generation of experiments 0, Q aims to 
trap H atoms as this is likely necessary for precision CPT 
and gravity tests. Neutral H atoms have a small per- 
manent magnetic moment, and can be trapped in the 
magnetic minimum of a so-called Minimum-B trap 0. 
The magnetic minimum can be created by two axially 
separated mirror coils which create an axial minimum, 
and a multipole field, such as an octupole |7|, l8|, which 
creates the radial minimum. In all current schemes, the 
Minimum-B and Penning-Malmberg traps must be co- 
located because the p's, e+'s, and Hs must all be trapped 
in the same spatial region. Thus, in cylindrical coordi- 
nates (r, ^, z), the net magnetic field will be 



f 008(46*) -^sin(46') +'Bu{r,e,z) 

(1) 



when using an octupole. Here Ryo is the trap wall radius, 
5^ is the octupole field at the wall, and Bm(^, z) is the 
field of the mirror coils. The mirror coils were not en- 
ergized for the data taken for this paper; henceforth we 
will set Bm = 0. 

Minimum-B traps are shallow (of order 0.7 K/T per 
Bohr magneton), and experimentalists have not yet 
learned to synthesize H with sufficiently low energy to 
be trapped. One obstacle to progress has been the lack 
of detailed information about the p cloud [10] dimensions. 
Until recently, only two techniques that measure the p ra- 
dial profile have been reported in detail. The first, based 
on p annihilation on the background gas [11], yields a 
crude [~ 4 mm (Icr)] three-dimensional image of the p 
cloud. To observe a sufficient number of annihilations, 
the background gas pressure must be much higher than 
is normally used when synthesizing antihydrogen atoms. 
This may influence the p cloud dimensions. The sec- 
ond interpolates the density profile from two destructive 
measurements [12]: the total p number, and the number 
that are located within a fixed radius set by an aperture. 
The reconstruction makes assumptions about the appli- 
cability of the global thermal equilibrium state of these 
plasmas [l3|, [m], and about the p temperature. We note 
that with our diagnostics (reported here and in [9]) we 
have seen many long-lived radial profiles that are not in 
global thermal equilibrium. 

Recently we described a diagnostic that gives high 
quality information about the radial profile. The diag- 
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FIG. 1: (Color online) Schematic diagram of the ALPHA apparatus. Particles are confined axially in an electrostatic well 
formed by biasing cryogenically- cooled, cylindrical electrodes centered on the trap axis. The axial magnetic field, graphed 
below the schematic, confines the particles radially. The p's were caught with the inner solenoid on, in a field of 3 T, as shown 
by the blue dashed curve. The inner solenoid was ramped off before transfer of the p's to the mixing region. The experiments 
described here were done in the 1 T field shown by the red solid curve. The MCP/Phosphor screen used to take images 
of the inner regions of the e~ plasmas and p clouds is located to the right of the parts of the apparatus shown here, in a field 
of 0.024 T. 



nostic is based on a MCP-phosphor screen system [9^ . (A 
similar system has also been reported by the ASACUSA 
collaboration [l5[.) Unfortunately, apertures limit the 
size of the p cloud that we can measure with our MCP- 
phosphor system; typically we cannot measure the profile 
beyond radii of 1.5-3.0 mm, depending on the local mag- 
netic field in which the p's are trapped. Some p clouds 
are completely imaged by this system, but others are far 
larger, and can extend all the way out to the walls of our 
trap at radius = 22.3 mm. Here, using the ALPHA 
collaboration trap |3] , we describe a new diagnostic that 
probes the outer radial profile based on measurements of 
ballistic [l^ losses induced by an octupole magnet. Af- 
ter a brief description of how we load particles into the 
trap, we describe the diagnostic. Then we discuss tests 
used to validate its performance, and close with several 
examples illustrating its use. 



degrading foil. The p's are cooled by collisions with a 
pre-existing electron (e~) plasma [17]. Multiple p pulses 
can be caught and cooled, each adding about 40,000 p's 
to the trap. Typically we use four such "stacks" in the 
data presented here. The e~ plasma is then ejected by 
fast manipulations of the electrostatic well that leave the 
massive p's behind. After cooling and e~ ejection, the 
p's are transferred, via manipulations of the electrostatic 
well potentials, to the "mixing" region of the trap. The 
octupole magnet [8] we use to determine the p radial pro- 
file is centered over this region. Positrons, when needed, 
are transferred from our Positron Accumulator p^, [loj 
and recaptured in the region indicated in Fig.[TJ They are 
then transferred to the mixing region via manipulations 
of the electrostatic well potentials. 



III. DIAGNOSTIC DESCRIPTION 



II. TRAP LOADING CYCLE 

We load our trap by accepting a pulse of p's from the 
AD. The p's enter the apparatus from the left (see Fig.[T]), 
and are slowed in a degrading foil. They reflect from a 
repelling potential at the far end of the "catching" region 
of the trap, and are then captured into an electrostatic 
well by quickly erecting an electrostatic barrier, at the 
near end of the trap, before they can escape back to the 



To understand how the radial diagnostic works, it is 
helpful to visualize the field lines from the solenoid and 
octupole coils. The field lines originating from a cir- 
cular locus of points in the plane transverse to z form 
four-fluted cylindrical surfaces; the flutes at each end are 
rotated by 45° with respect to each other. An exam- 
ple of the resulting surfaces is shown in Fig. [2l Fig. [3] 
shows an image of one quadrant of the field lines, gener- 
ated by passing e~'s through the octupole and onto our 
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FIG. 2: (Color online) Magnetic field from the octupole and 
solenoid coils. The vectors on the left represent the directions 
of the axially- invariant field from these coils. The surface is 
created by following the field lines from a radially centered 
circular locus; the lines shown within the surface are field 
lines. 
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MCP/Phosphor screen 0. 

Antiprotons confined by the electrostatic well within 
the octupole bounce back and forth while following the 
magnetic field lines |20|. Antiprotons that are on field 
lines that extend to the physical trap wall before reach- 
ing one of the electrostatic walls will follow them there 
and annihilate. For a given end-to-end bounce length 
L, field lines lying outside of a critical radius Vc at the 
trap center will hit the wall, while those lying inside the 
critical radius will not. The normalized critical radius is 



This relation is depicted in Fig. [H The longer the trap, 
and the stronger the octupole field, the smaller the crit- 
ical radius. The normalized critical radius is never very 
small because the octupole field, which scales as r^/i?^, 
is very weak near the trap axis relative to its strejigth at 
the wall. This is advantageous for confinement [7], as a 
large cloud survives and the inner core of the p cloud is 
not strongly perturbed by the multipole field. However, 



Octupole Off Octupole On 




FIG. 3: (Color online) Field lines imaged by passing a cir- 
cular e~ plasma through the octupole with the octupole off 
and on. Apertures [9] form the image boundaries and limit 
us to viewing only one quadrant of the octupole field map. 
The distortion evident in the right-hand image corresponds 
to one of the flutes at the end of the magnetic surface shown 
in Fig. m 
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FIG. 4: (Color online) The normalized critical radius [Eq. 
as a function of the octupole strength Bw and orbit length L. 
The alternate axes shown at the top isolate the dependence 
on each parameter while holding the other fixed at a typical 
value. 



as we show below, it limits the observable minimum p 
radius to about 7 mm for a 135 mm long well. If we had 
used a quadrupole instead of an octupole, we could have 
measured radial distributions to much smaller radii; for 
instance, to 0.24 mm for equivalent parameters. Such a 
small critical radius would be very useful as a diagnostic, 
but could make it difficult to synthesize H. 

The ballistic loss of particles on trap walls in the pres- 
ence of a multipole field was first identified with elec- 
trons in a quadrupole magnet [l^ . This process is easier 
to study with p's than with e~'s, however, because indi- 
vidual p annihilations can be detected and localized on 
the trap wall with a position sensitive detector. The de- 
tector j23| comprises three layers of silicon cylindrically 
arrayed around the trap axis just outside of the octupole 
magnet (see Fig. [1]). It is not yet fully deployed, but, 
using a partial system consisting of 10% of the full sys- 
tem, we observe (Fig. [5j) that p's hit the wall at the ends 
of the electrostatic well. We expect to observe this type 
of loss pattern as it is at the ends of the trap that the 
accessible field lines extend furthest outward; we note, 
however, that annihilations tend to occur at the ends of 
the electrostatic well even in the absence of an octupole 
field [11]. 

For the experiments reported in Fig. [6 HT31 annihila- 
tions were detected by scintillators coupled to Avalanche 
Photo Diodes (APDs). As with the silicon detector, the 
scintillators are cylindrically arrayed around the trap axis 
just outside of the octupole magnet. Annihilations are 
identified by the firing of more than one scintillator in a 
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FIG. 5: (Color online) Axial positions where the p's hit the 
trap wall under the influence of the octupole. The horizontal 
bar indicates the the axial extent and position of the electro- 
static well confining the p's. The loss is greatest near the ends 
of the confining electrostatic well. The positions are deter- 
mined by a position-sensitive particle detector which monitors 
the p annihilation products; the vertical lines at z = ±115 mm 
indicate the axial extent and position of the detector. 



150 ns coincidence window, and we detect annihilations 
with greater than 50% efficiency. The detector back- 
ground noise is of order a few events per second. Timing 
modules correlate annihilations with experimental opera- 
tions and conditions such as the strength of the octupole 
field. 

To measure the size of a p cloud, we first transfer it 
into an electrostatic well in the octupole field region; the 
octupole field is turned off during the transfer. We then 
measure the p kinetic energy by monitoring the rate at 
which the p's escape as we slowly lower one endwall of 
the electrostatic well [24] . Typically we find that the en- 
ergy is between 1 and 15 eV; the energy depends on the 
details of the transfer process and the electrostatic well 
potentials. This measurement is destructive, but since 
the energy is largely set by the electrostatics, not by the 
p radial profile, it is sufficient to measure this energy 
once for a series of profile measurements. From this en- 
ergy, we determine the bounce length L of the p's in the 
electrostatic well. The uncertainty (and spread) of the p 
energy sets the uncertainty in the orbit lengths quoted 
in the figure captions. Finally, for each p cloud that we 
want to analyze, we slowly ramp up the octupole field 
Byj while monitoring the losses. From the time history 
of the losses, we can invert Eq. ([2]) to reconstruct the 
radial distribution of p's: 



N{t) 



(3) 



Here B^{t) is the octupole field at time t, rc[Bw{t)] is the 
instantaneous critical radius, and dVc/dB^ is evaluated 
at the instantaneous field B..,. The raw data from our 



detector is binned in intervals of time Ato = 1 ms; we re- 
bin the data into intervals ranging between At = 0.333 s 
(45 s and shorter octupole ramp times) and 1.332 s (180 s 
ramp times) to decrease the scatter. N{t) is the number 
of counts in the bin centered around t. The mapping de- 
fined by Eqs. ([2]) and (|3]) is nonlinear; points are closer 
together in r at small radii than at large. To further 
reduce the scatter at small r we rebin n(r) so that the 
spacing between successive points in r is never less than 
0.075 mm. 



IV. VALIDATION TESTS 

Typical data are displayed in Fig. [6l which shows 
the radial profile of two otherwise identically prepared 
p clouds stored in wells of different length. Changing 
the well length should not change the radial profile of 
identically prepared p clouds, and as expected, the mea- 
sured profiles are almost identical over their common 
range. However, as predicted by Eq. ([2]), changing the 
well length does change the minimum radius observable 
with the diagnostic from about 7.0 mm for the 135 mm 
well, to 9.6 mm for the 65 mm well. 

Figure [3 compares the radial profiles of identically pre- 
pared p clouds held in a flat-bottomed well, and in a 
nested well similar to those used to synthesize H [1]. The 
well length inferred from the measured p energies was 
130 mm for the nested well, which is slightly shorter than 
the 135 mm length inferred for the flat well. Changing the 
well shape should not change the radial profile because 
the azimuthally-symmetric electrostatic well fields do not 
induce radial transport. As expected, the measured pro- 
files are nearly identical. Thus, the diagnostic is indeed 
independent of the well shape so long as the proper well 
length is employed in the analysis. 

As the octupole ramps, outward diffusion m- 
creases for those p's that are still within the critical ra- 
dius; if this diffusion were too fast, the profiles would be 
suspect. We have established that the diffusion is not 
fast on the time scale of the octupole ramp by compar- 
ing (Fig. [8|) the radial profiles of identically prepared p 
clouds taken with ramps of 45 (our standard ramp), 90, 
and 180s. The differences between the curves are not 
large. 

The diagnostic described here would have little utility 
if all reconstructed radial profiles were identical; Figured 
shows that radial profiles of p clouds that are differently 
prepared can be dissimilar. Figure [9] also shows that the 
load-to-load reproducibility of the p profiles is quite good. 

Measurements taken with our MCP/phosphor screen 
diagnostic confirm that the central density is not sig- 
nificantly perturbed by cycling the octupole field. For 
instance, for parameters identical to the nested well pro- 
file shown in Fig. [71 the total number of p's within the 
MCP/phosphor apertures varied by less than 4% on two 
successive shots, one with the octupole off and one with 
it ramped up and then back down. This discrepancy is 
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FIG. 6: (Color online) Comparison of the radial profiles of 
otherwise identical p clouds held in wells of different length. 
Panel a) shows the electrostatic well potentials ^{z) for the 
two cases; the horizontal bars indicate the axial extent and 
position of the p orbits before the application of the octupole 
field. Panel b) shows the time history of the p annihilations as 
the octupole field is ramped up. Panel c) shows the resulting 
radial profiles. In all graphs, the green solid curve corresponds 
to the longer well (135 ± 5 mm) and the red dash curve cor- 
responds to the shorter well (65 =b 5 mm). The maximum Bw 
at the end of the 45 s ramp was 1.54 T, and Bz = 1.03 T. 
At the inner radii, Eq. (|2]) predicts that the ±5 mm length 
uncertainty/ spread engenders a radial uncertainty of about 
±0.12 mm at 135 mm, and ±0.30 mm at 65 mm. Near the 
wall, the uncertainty predicted by Eq. (|2| diminishes, but the 
time binning engenders an uncertainty of about ±0.25 mm. 
The error bars indicate the size of the typical calculated sta- 
tistical error. Both p clouds were collected with four stacks. 



well within the shot-to-shot variation of our loads. This 
result, taken together with the results shown in Figs. [6]- 
m establish that ramping the octupole field is a robust 
method of obtaining the radial profile that is largely inde- 
pendent of the details of the ramp speed and well shape. 



V. OBSERVATIONS 

We have used our new diagnostic to characterize our p 
manipulation sequences, and to study interesting physics 
issues. In this section, we outline four of these measure- 
ments; all need further study. 

As described earlier, we can stack multiple p pulses 
from the AD. Figure [10] shows the p profile for two, 
three, and four stacks. The stacks add to each other 
without significantly changing the radial profile. The re- 
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FIG. 7: (Color online) Comparison of the radial profiles ob- 
tained with flat (green solid) and nested well potentials (red 
dash). The well lengths were 135 ± 5 and 130 ± 5 mm respec- 
tively. The graph descriptions and all other parameters are 
the same as in Fig. [6] 
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FIG. 8: (Color online) Comparison of the radial profiles with 
octupole ramps of 45 (green solid), 90 (red short-dash), and 
180 s (blue long-dash). The well length in each case was 130 ± 
5 mm. The graph descriptions and all other parameters are 
the same as in Fig. [6] 
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FIG. 9: (Color online) Radial profiles for two sets of p clouds 
that were prepared differently; the e~ cooling plasmas used for 
the two sets came from different e~ sources. The figure also 
shows that the load-to-load reproducibility of the p clouds 
is high; the set labeled I compares two loads, while the set 
labeled II compares three. The AD and our apparatus can be 
quite reproducible; the two profiles in set I were measured 23 
hours apart on different AD shifts. (Note that the clouds were 
analyzed in different shape wells, one (I) of length 135 it 5 mm 
in a flat well, and the other (II) of length 130 it 5 mm in a 
nested well. The ramp time for set II was slightly shorter than 
for set I: 36 s instead of 45 s. However, as verified in Figs. [71 
and [H these difference should not affect the radial analysis. 
Finally, only two stacks were used in set II; the profiles for this 
set were normalized to four stacks.) The graph descriptions 
and all other parameters are the same as in Fig. [6l 



suits obtained when only one stack is accumulated are 
quite different, however. The profile is completely con- 
tained within a radius of 7 mm and is not visible with 
this diagnostic. We suspect that the difference is due to 
straggler e~'s from the degrader accidentally captured 
during the first (and subsequent) p injections. These 
e~ 's are captured by the same electrostatic well manipu- 
lations used to capture the p's. After capture, they cool 
and thermalize via cyclotron radiation and collisions, and 
join the deliberately captured cooling e~ plasma; we ob- 
serve that the number of e~'s in this plasma increases 
with the number of stacks. The straggler e~'s are likely 
emitted from the degrader over the entire area hit by the 
p's, and, if the radius of this area is greater than the ra- 
dius of the deliberately injected e~ plasma, the plasma 
radius will increase. This will increase the size of the 
captured p cloud [9] . It will also increase the fraction of 
the degraded p's captured [9]; we observe this fraction 
increasing from about 45% on the first stack to over 90% 
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FIG. 10: (Color online) Comparison of the radial profiles for 
two (blue long-dash), three (red short-dash), and four (green 
solid) stacks. The well length was 130 =b 5 mm. The graph 
descriptions and all other parameters are the same as in Fig.[6l 



on later stacks. 

The transfer process from the catching region of our 
trap to the mixing region leaves the p's situated in a 
short well on one side of the final trapping well. From 
this short well, the p's are injected into the final well. 
Normally, we do this gradually, by smoothly changing the 
potentials over a 1 ms time period. When we change the 
potentials abruptly, on a time scale of approximately 3 //s, 
p's are lost on injection, and the p cloud's radius increases 
significantly, as shown in Fig. [Til There is no obvious 
mechanism for the immediate loss and cloud expansion. 

Figure [12] shows the very different radial profile ob- 
tained when we do not eject the e~'s before transfer and 
analysis. The antiprotons form a hollow ring around 
the trap center. This type of distribution is compati- 
ble with the global thermal equilibrium of a mixed e~-p 
plasma, which places the p's in a halo surrounding the 
e~ plasma [26] when the particles are sufficiently cold. 
However, we observed losses during the transfer process 
that could have preferentially hollowed the distribution 
and produced the observed profile. 

Note that the p's likely cool via collisions with the e~ 's 
during the octupole ramp. This would shorten the axial 
extent of the p orbits, and thus introduce some uncer- 
tainty into the reconstruction of the radial profiles via 
Eq. (|3|) as it introduces variation in L. This is particu- 
larly true if the p's cool into the side wells, where their 
orbit length would decrease abruptly by more than a fac- 
tor of two. This effect would cause us to erroneously 
reconstruct, via Eq. (|3|), some charge to be at falsely low 



7 



2x103 




-100 100 

z (mm) 



20 40 

Time (s) 



O 




10 15 
r (mm) 



FIG. 11: (Color online) Comparison of the radial profiles ob- 
tained with a gentle (green, solid) and abrupt (red short-dash) 
injection into a long well. The blue dash curve in Panel a) 
shows the pre-injection well structure (the p's start in the 
leftmost well) and the green solid curve shows the final well, 
which has a length 135 ± 5 mm. The graph descriptions and 
all other parameters are the same as in Fig. (6] 



radii, probably below the 7 mm radius visible to us with 
this diagnostic. Thus, cooling does not explain the halo 
visible in Fig. [121 This very interesting result needs fur- 
ther study. 

Finally, in Fig. [131 we show radial profiles for a mixed 
e+-p plasma. As the density of the e+ plasma is in- 
creased, p's appear to be transported outward. Here, as 
described in the previous paragraph, the interpretation 
of the results is complicated by cooling of the p's (on the 
e+ in this case.) Cooling will again cause some charge to 
appear at falsely low radii, and this very likely causes us 
to underestimate the outward movement of the p's. 

A possible explanation of the outward movement 
shown in Fig. [131 is that it is the result of the forma- 
tion of highly excited H that is either 1) ionized at the 
radial edge of the e+ plasma by its self consistent electric 
field, which is strongest at the edge, or 2) ionized by the 
vacuum electrostatic well fields. Note that the p's from H 
that was ionized within the e+ plasma radius would have 
the opportunity to recombine into H again, while those 
at larger radii would orbit unperturbed. With time, the 
p's remaining in the e+ plasma would be swept out to 



larger radii. Unpublished simulations of realistic antihy- 
drogen formation/field ionization cycles, using the code 
described in [l^, found similar transport. We do not yet 
have any other direct experimental evidence that this cy- 
cling is occurring. 
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FIG. 12: (Color online) Comparison of the radial profiles with 
and without electrons (dashed red and solid green lines, re- 
spectively. The well length was 130 di 5 mm, and the ramp 
time was 180 it 5 s. The graph descriptions and all other pa- 
rameters are the same as in Fig. [6l 



VI. CONCLUSIONS 



We have shown that we can determine the outer ra- 
dial profile of p's stored in a Penning-Malmberg trap by 
monitoring the losses induced by ramping an octupole 
magnet. This technique complements direct imaging of 
the inner radial profile , and provides more precise and 
reliable information than earlier techniques [ijj, [l2| . We 
have tested the diagnostic by varying the electrostatic 
well length and shape, and by varying the ramp time, 
and we have used the diagnostic to study several pro- 
cedures and manipulations pertinent to the synthesis of 
antihydrogen atoms. 

This work was supported by CNPq, FINEP (Brazil), 
ISF (Israel), MEXT (Japan), FNU (Denmark), NSERC, 
NRC/TRIUMF (Canada), DOE (USA), EPSRC and the 
Leverhulme Trust (UK) and HELEN/ ALFA-EC. 
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FIG. 13: (Color online) Comparison of the p radial profile 
with different density positron plasmas. The green, solid 
curve shows the profile with no e^, and the red short-dash 
and blue long-dash curves show the profile with 13 million and 
25 million e^ respectively. The well length was 85it 5 mm, the 
maximum field was 1.20 T, and the ramp time was 20 s. Only 
one stack was captured, but the e~ cooling plasma was cre- 
ated by a secondary e~ source which makes a large radius 
plasma; thus, unlike in Fig. [TOl p's are visible with only one 
stack. The graph descriptions and all other parameters are 
the same as in Fig. [6] 
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